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We have carried out a detailed experimental investigation of the static properties of planar Joseph- 
son tunnel junctions in presence of a uniform external magnetic field applied in an arbitrary orien- 
tation with respect to the barrier plane. We considered annular junctions, as well as rectangular 
junctions (having both overlap and cross-type geometries) with different barrier aspect ratios. It 
is shown how most of the experimental findings in an oblique field can be reproduced invoking the 
superposition principle to combine the classical behavior of electrically small junctions in an in-plane 
field together with the small junction behavior in a transverse field that we recently published [R. 
Monaco, et al,, J. Appl. Phys. 104, 023906 (2008)]. We explore the implications of these results 
in supposing systematic errors in previous experiments and in proposing new possible applications. 
We show that the presence of a transverse field may have important consequences, which could be 
either voluntarily exploited in applications or present an unwanted perturbation. 

PACS numbers: 74.50. +r 



I. INTRODUCTION 



One of the earliest experiments involving Josephson junctions and magnetic fields has been the measurement of the 
magnetic diffraction pattern^ i.e., the dependence of the junction critical current I c on the amplitude of an externally 
applied magnetic field H a . Traditionally, since the discovery of the Josephson effect in 1962, the magnetic diffraction 
pattern I c (H a ) of planar Josephson tunnel junctions (JTJs) has been recorded with the magnetic induction field 
applied in the junction plane to avoid the huge computational complications of taking demagnetization effects into 
account, when a transverse magnetic component is present. A number of important results have been derived from 
experiments under these assumptions - a prominent example being the determination of the London penetration depth 2 
\l from which one derives the Josephson penetration deptfP A j which sets the JTJ electric length scale. Nowadays, 
every textbook on the Josephson effect deserves at least one chapter to the magnetic diffraction phenomena. The 
simplest case is that sketched in Fig [I] of a rectangular JTJ placed in a uniform and constant external magnetic field 
parallel to one of the barrier edges. Let us choose the coordinate system such that the tunnel barrier lies in the z = 
plane and let 2L and 2W be the junction dimensions along the x and y- direct ions, respectively. Finally, let us assume 
that the JTJ is electrically small, meaning that its dimensions are both smaller than the Josephson penetration depth 
(2L,2W < A j) (absence of self-fields) and that its Josephson current density Jj is constant over the barrier area. 

If the externally applied field H a is along the x direction, H a = H x x, then the magnetic field H inside the junction 
is constant and equal to the external value, i.e., H = (i^,0,0). By integrating the Josephson equatiorP relating the 
Josephson phase <p to the magnetic induction field in the barrier H: 



Vcf>(x,y) = — ^H(^)xz, (1) 

in which d e is the junction magnetic thickness, fio is the vacuum permeability and <E>o = h/2e is the magnetic flux 
quantum, we readily obtain the spatial dependence of the Josephson phase: 

<j> = tvH x y, (2) 
with ft = 2ird e fio/<&o. Eq.Q leads to the well known Fraunhofer-like magnetic diffraction pattern 3 : 



I C (H X ) — Iq 
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where Iq = AJjWL is the zero field junction critical current and H c = $>o/2[iod e L is the so-called (first) critical field, 
i.e. the smallest field value for which the Josephson current vanishes. Barone and Paterno 4 generalized Eq.([3| to the 
case of an arbitrary orientation of the external magnetic field in the junction plane, H a = H x x + H y y. In such a case, 
still H = H a and the resulting magnetic diffraction pattern will be: 



Ic(Hx, H y ) — Iq 
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with H cx = ^o/2jaod e L and H cy = <&Q/2jiod e W . Unfortunately, the last equation cannot be easily generalized to the 
case of an arbitrary applied field orientation H a = H x x + H y y + H z z, simply because, when H z ^ 0, then H 7^ H a . 
The effect of a transverse magnetic field has been first considered in 1975 by Hebard and FultorPin order to provide 
a correct interpretation to some experimental data published in the same yeaiP. They observed that a transverse 
applied field H a = H z z induces Meissner surface demagnetizing currents j s feeding the interior of the junction and so 
generating a magnetic field H in the barrier plane such that j s = z x H. The problem of finding the j s distribution in 
a single superconducting field subjected to a transverse magnetic field has been analytically solved for different film 
geometries^. However, a planar JTJ is made by two overlapping superconducting films separated by a thin dielectric 
layer and the Meissner current distributions on the interior surfaces (top surface of the bottom film and bottom 
surface of the top film) require numerical approaches even for the more tractable electrode configurations. Recently 9 , 
the magnetic field distribution H x in the barrier of small planar JTJs has been numerically obtained in the case when 
an external magnetic field is applied perpendicular to the barrier plane, H a = (0,0, H z ). The simulations allowed 
for heuristic analytical approximations for the Josephson static phase profile (j) 1 - from which the dependence of the 
maximum Josephson current I C (H Z ) on the applied field amplitude was calculated for the most common electrode 
geometrical configurations (overlap, cross and annular junctions). Unfortunately, the theoretical findings could not 
be tested against experimental results due to the insufficiency of data available in the literature. 

One of the aims of this paper is to fill this vacancy. We have measured the transverse magnetic diffraction patterns 
of several planar JTJs with the most common geometrical configurations and compared the results with their expected 
counterparts. More generally, we have recorded the I c (H a ) when the applied field is oblique, that is, has non-zero in- 
plane and transverse components. To avoid complications and without loss of generality, we have chosen the in-plane 
component to be along one of the electrode axis, more specifically, along the x-direction: 



H a H x x + H z z, 

so that, as shown in Fig{l] the applied field forms an angle a with respect to the x-y plane, that is: 



(5) 



H x = H a cos a and H z = H a sin a, 



(6) 
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TABLE I: Relevant electrical (at T=4.2K) and geometrical parameters of the rectangular and annular Nb/Al ox /Nb Josephson 
tunnel junctions quoted in this paper. The Josephson critical current density was J c = 3.9 kA/cm 2 (corresponding to Aj ~ 6 fim) 
for all samples, except for sample #C having J c = 80 A /cm 2 (Aj c± 42 fim). The experimental results obtained for these samples 
will be presented in Sec |III| as follow: overlap-type junctions in Sec |III A| cross-type junctions in Sec |IIIB| and annular junctions 
in Sec HlTCl 



with H a = yj H 2 + H 2 . We will demonstrate that the experimental oblique magnetic diffraction patterns can be 
nicely reproduced by properly extending the theoretical framework of Ref. 9 . This paper is constructed as follows. In 
Secjn|we will present the samples used for the measurements and describe the experimental setup. Sec |III| will report 
on the exper imental results obtained for tho se sam ples whose barrier has a rectangular shape (overlap-type junctions 
in Sec |III A and cross- type junctions in Sec III B). Sec III C will be devoted to the annular JTJs. Then, in Sec IV 



we will discuss how to generalize the theoretical analysis of the effect of a transverse magnetic field to the case of 
an oblique field. Finally, the discussion and the interpretation of the measurements will be given in Sec(V| while the 



conclusions will be presented in Sec VI 



II. THE SAMPLES 



High quality Nb/Al — Al ox /Nb JTJs were fabricated on 0.35 mm thick silicon substrates using the trilayer technique 
in which the junction is realized in the window opened in a S1O2 insulator layer - details of the fabrication process 
can be found in RefP^l The so called passive or idle region, i.e. the distance of the barrier borders to the electrode 
borders, was on the order of 1-2 \im for all the junctions. The thickness of the S1O2 insulator layer was 400 nm. The 
demagnetization currents strongly depend on the electrode thicknesses relative to the London penetration depth. For 
our samples the nominal thicknesses of the base, top and wiring Nb layers were 200, 100 and 500 nm, respectively. 
Considering that the London penetration depth for Nb film is A^ ~ 90 nm 2 , we see that our samples satisfy the 
thick film approximation. For all samples the high quality has been inferred by a measure of the I-V characteristic at 
T = 4.2 K. In fact, the subgap current I sg at 2mV was small compared to the current rise AI g in the quasiparticle 
current at the gap voltage V g , typically AI g > 20I sg ; the gap voltage was as large as V g = 2.8 mV. The geometrical 
and electrical (at 4.2 K) parameters of the seven samples quoted in this paper are listed in Table I. For the rectangular 
junctions #A-F, beside their dimensions 2L and 2W along the x and ^-directions, respectively, we also report the 
junction aspect ratio (3 = L/W. (As shown in RefP, this geometrical parameter turns out to be crucial for the magnetic 
field line distribution in the barrier of a JTJ subjected to a transverse magnetic field.) All the samples belonged to 
the same fabrication batch (except sample #C). Let us observe that for the overlap- type junctions -j^A and j^B the 
zero field critical current Io was as large as the theoretical value 0.7 AI g predicted for strong-coupling Nb-Nb JTJs, 
indicating the absence of self- field effects. The critical current density has been calculated as 11 J c = 0.7 AI g /A in which 
AI g is the measured quasiparticle current step at the gap voltage and A is the junction nominal area (A = 2L x 2W 
for rectangular junctions and A = 7r(r 2 — rf) for the annular junction). The Josephson critical current density was 
J c = 3.9/c A/ cm 2 for all samples, except for sample #C having J c = SO A/ cm 2 . The values of the barrier magnetic 
thickness d e = 2A^ ~ 180 nm has been used to calculate the Josephson penetration depth Aj = ^/ V>o/27r/io^e^c- (I n 
the thin film limit, Aj can be better determined by using the expression for d e found by Weihnacht 12 .) Accordingly, 
all samples had A j c± 6 /xra, except sample #C which had A j c± 42 \im. In other words, as far as the electrical length 
concerns, all samples can be classified as intermediate length junctions (2L,2W ~ Aj), except sample #C that is a 
long (2L » Aj) unidimensional (2W « Xj) overlap-type JTJ. 

We now come to the definition of the parameters A^, and their ratio tjr whose experimental values are 
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FIG. 2: (Color online) Definition of the parameter Ar as the width of the magnetic field range in which I c (H a ) > (2/ti)Iq. By 
definition, for a Fraunhofer-like magnetic diffraction pattern, Ar coincides with the junction (first) critical field H c . 



reported in Table I for the rectangular junctions (A^, A^ and tja for the annular junction). As already mentioned 
in the Introduction, it's well known that the magnetic diffraction pattern of an electrically small rectangular JTJ in 
the presence of an in-plane field perpendicular to one of the barrier edge follows the Fraunhofer pattern in Eq.([3| 
characterized by a periodic amplitude modulation. As depicted in Fig(2j the value H c of the applied field where first 
the critical current vanishes is called the (first) critical field. (H c is a measure of the response (for small field) of 
the junction critical current to the applied field: the smaller is H C1 the larger is the junction response.) For those 
samples whose I c (H a ) is still amplitude modulated, but follows a different pattern (for example, annular, circular 
and rhombic junctions), the critical field H c can still be defined as that value of external field H a where first the 
critical current nulls, I C (H C ) = 0. Further, for those samples whose I c (H a ) shows modulation lobes, but never 
vanishes (as examples, small junctions with nonuniform tunneling current and long JTJs), the critical field can still 
be obtained extrapolating to zero the first modulation lobe. However, in those cases in which the critical current I c 
is a monotonically decreasing function of the applied field H ai the concept of critical field H c looses its meaning and 
a new feature has to be introduced to characterize the behavior of I c (H a ) for small fields. A theoretical example is 
offered by a Gaussian shaped junction subjected to an in-plane magnetic field that is characterized by a Gaussian 
magnetic diffraction pattern^. A practical example is given by a square cross junction in a transverse field, whose 
I c (H a ) decreases with H a with no measurable modulation 14 . The experimental magnetic diffraction patterns that will 
be reported in the next section span all kinds of behaviors from Fraunhofer- type to 1/iJ^-like (with v > 0). Therefore, 
as the new and universal figure of merit to characterize the response of the critical current to the externally applied 
field amplitude, we have chosen the width of the magnetic field range Ar in which I c (H a ) > (2/ty)Iq c± 0.64/ 
(see Fig {2]). Considering that, when Eq.Q holds, I c (H c /2) = (2/tt)Iq, the value of the prefactor stems from the 
requirement that the new merit figure Ar numerically equals the critical field H c whenever the measured magnetic 
pattern follows a Fraunhofer dependence, i.e., Ar = H c . In all other cases, generally speaking, Ar 7^ H c . In our 
notation, A^ and A^ are the merit figures of, respectively, an in-plane (a = 0) and transverse (a = 90°) magnetic 
diffraction pattern. With a similar reasoning, we define the parameter A^ for annular junctions as the width of 
the magnetic field range Ar in which I c (H a ) > /iio, with \i ~ 0.67. The slightly different prefactor stems from the 
fact that the in-plane diffraction pattern of a small annular junction (with no trapped fluxon) follows a Bessel-type 
dependence 15 : I c (H a ) = Io \Jo(CiH a /H c )\ in which Jo is the zero-order Bessel function and £i w 2.405 is its first zero. 
Now H c = $0/ )iiod e C, where C is the ring mean circumference. 

The measurement of A requires an external field smaller than the one required for H c ; henceforth, this new parameter 
also turns out to be a very useful quantity whenever the junction critical field cannot be experimentally determined 
since it exceeds the irreversible field, i.e., when the Abrikosov vortices first enter into the superconducting films and 
become pinned in the junction^. For our samples the transverse irreversible field was about 5mT (50 Gauss). 
The ratios tjr = A^/A^ and tja = A^/A^ provide a direct comparison between the I c response to a transverse field 
relative to the in-plane field; specifically, 77 < 1 means that the junction critical current modulates faster when the 
applied field is transverse. In a recent paper^we already provided an experimental proof that a transverse magnetic 
field can be much more capable than an in-plane one to modulate the critical current I c of a planar JTJ with proper 
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barrier and electrodes geometry requirements. This property was first obtained and exploited in the context of a 
detailed investigation of the phase symmetry breaking during fast normal-to-superconducting phase transitions of 
long annular JTJs^. 

Our setup consisted of a cryoprobe inserted vertically in a commercial LHe dewar. The cryoprobe was magnetically 
shielded by means of two concentric magnetic shields: the inner one made of Pb and the outer one of cryoperm. Inside 
the vacuum tight can of the cryoprobe, a non-magnetic insert holds a chip mount with spring contacts to a Si chip 
with planar JTJs. With reference to the coordinate system in Fig{l] the chip was positioned in the center of a long 
superconducting cylindrical solenoid whose axis was along the x-direction (within less than 1° of accuracy) to provide 
an in-plane magnetic field. In order to provide a transverse magnetic field, a superconducting cylindrical coil was 
placed 5 mm far from the chip with its axis oriented along the ^-direction (within less than 3° of accuracy) . Two 
independent low-noise dc current sources were used to feed the solenoid and the coil in order to expose our samples 
at magnetic fields having arbitrary magnitude and orientation (in the a>z-plane). The field-to-current ratio was 
3.9 iiT/mA for the solenoid and 4.4 fiT/mA for the coil. These values have been numerically obtained from Comsol 
Multiphysics magnetostatic simulations in order to take into account the strong correction to the free-space solution 
due to the presence of the close fitting superconducting shield^. 




FIG. 3: (Color online) Tridimensional plot showing for the square overlap junction [[3—1) the recorded magnetic diffraction 
patterns I c (H a ) for different values of the field orientation a (with Aa = 10°). 



III. THE MEASUREMENTS 

In this Section we present the experimental oblique magnetic diffraction patterns relative to planar JTJs having 
the seven different electrode configurations listed in Table I. Subsection A, B and C will concern samples having, 
respectively, overlap, inline and annular geometry. The theoretical interpretation of our data-sets will be given in the 
next Section. The angle a that the external oblique field forms with barrier plane could be experimentally spanned 
in the interval [— 7r,7r]; however, according to Eq.(|6|, an angle rotation of ±7r is equivalent to an inversion of the 
field direction, i.e., of the field amplitude H a — > —H a . For this reason we will only present data for a in the [0,7r] 
interval with the amplitude H a assuming both negative and positive values. Further, by denoting with /+ and I~ 
the positive and negative critical currents, respectively, we always had I~{H a ) = /+(— if a ), as expected, due to the 
absence of any measurable stray fields in our setup. For this reason, we will only present data for /+, which we will 
simply call I c . We stress that, in recording the I c vs. H a curves, we took special care that the applied field never 
exceeded the reversible field, so it was not expected that the applied field penetrated the films. Furthermore, through 
measurements of the sample's I-V characteristic, it was verified that H a was so small as not to affect the energy gap. 
Finally, the raw experimental data were postprocessed to take into account the difference in the solenoid and coil 
field-to-current factors. 
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FIG. 4: (Color online) Magnetic diffraction patterns of the square overlap junction (/3 — 1) for different a values: a) a = 
(in-plane field), b) a = 50°, c) a = 90° (transverse field), and d) a = 140°. The experimental data are presented by closed 
circles. The solid lines, when present, are the best Fraunhofer fit, while the dotted lines are the results of the calculations 
described in Sec |IV| For a = 0, by construction, the calculations reproduce the Fraunhofer shape. 



A. Overlap-type junctions 

We begin with an intermediate length, square overlap-type JTJ, namely sample j^A in Table I (2L = 2W c± 1.6Aj). 
Figgis a tridimensional plot of the magnetic diffraction patterns recorded for different a values (with Aa = 10°). 
Since, for this sample, I c (—H a ) = I c (H a ), we only show the data for H a > 0. It is evident that I c (H a ) smoothly, 
but drastically, changes with the field orientation a. To be clearer, in Figs j4] we report the magnetic patterns for four 
selected a values. For a = the applied field is in the barrier plane (z = 0) and, as seen in Figj4^i, I c (H a ) closely 
follows a Fraunhofer- like behavior, as expected. The small discrepancy between the experimental data (closed circles) 
and the Fraunhofer fit (solid line) can be ascribed to the fact that junction dimensions are slightly larger than the 
Josephson penetration depth. Increasing a, in the beginning the junction critical field H c (or equivalently the width 
of the pattern main lobe Ar defined earlier) first slowly decreases until it reaches an absolute minimum when a 50° 
(see Fig j4]a), and later on quickly increases until it reaches an absolute maximum when a ~ 140° (see Figj4]i.) In 
FigjZJ: we also report the transverse (a = 90°) magnetic pattern to evidence how much it differs from a Fraunhofer 
dependence. It is worth stressing that, whenever a ^ Irrm (with integer m), the magnetic diffraction pattern looses 
the modulation periodicity H cn = nH c \ featuring the Fraunhofer behavior; more specifically, the distance between 
two adjacent minima increases as we move to larger fields, i.e., H cn > nH c i. Each plot in Figs J4] explicitly reports the 
corresponding position of the measured A#. 

The a dependence of A# (normalized to A^) for the sample j^A is summarized in Figj^i (solid circles). Ar(o) 
is reported in Figs|5}>d for the samples #5, #C, and #D, overlap-type JTJs having aspect ratios, respectively, 
(3 = 0.25, 0.08, and 4. The insets in Figsj5^i-d sketch for each sample its electrode configuration and its orientation 
with respect to the Cartesian coordinates chosen in the Introduction (see FigfTl). We used a vertical log-scale for those 
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samples having (3 < 1. Each plot in Figs|5]is characterized by an absolute maximum achieved when a = c^m- The 
olm values, quoted in the last column of Table I, were found to monotonically depend on the tjr ratios which, in turn, 
scale with the (3 ratios. In SecjVjwe will discuss a simple theoretical approach aimed to find the a dependence of A# 
and the relationship between «m and m?, as well. 

We like to point out that the only measurements similar to those reported in Figj5]can be found in a pioneering 
paper dated 1975 by Rosenstein and CherP. They measured the first and second junction critical fields in an oblique 
magnetic field for an overlap-type planar JTJ having (3 ~ 0.5 (and formed by two 300 nm thick Pb electrodes of 
unequal widths). They found that both H c \ and H C 2 reach their maximum values when the field orientation is about 
8° off the in-plane direction (% ~ 172° in our notation). This value is consistent with our findings. 
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FIG. 5: (Color online) Magnetic field range Ar vs. a (in degrees) for overlap junctions with different aspect ratios /3: a) (3 = 1, 
b) /3 = 0.25, c) /3 = 0.08, and d) (3 = 4. The experimental data are presented by closed circles. The open stars result from the 
calculations described in Sec |IV| while the solid lines, when present, are the result of a simple theory developed in Sec(V| The 
insets sketch for each sample its electrode configuration and its orientation with respect to the chosen Cartesian coordinates. 



B. Cross-type junctions 



Figure |6] shows the transverse magnetic pattern of the square cross junction (sample j^E in Table I) on a log- log 
plot. The inset displays the same data on linear scales. We observe that the critical current monotonically decreases 
as the external field is increased. The experimental data indicate that for large fields, I c oc H~ 2 , in contrast with the 
simple inverse proportionality suggested by Miller et aZH 

In Figs (7^) and b) we report the oblique magnetic diffraction patterns of the two cross- type junctions quoted in Table 
I, respectively, JJ#E and JJ#F. For these samples the in-plane patterns are skewed due to the self-field effects. 
The skewness is more pronounced for the asymmetric sample, having 2L « 3Aj and Iq = I c (H a = 0) = 3.6mA (of 
course Io does not depend on a). As we move from an in-plane field (say a = 0) to a transverse one (say a = 90°), 
the skewness gradually disappears and, keeping increasing a toward 180°, the skewness changes its polarity; in other 
words, J c (if a ,7r/2 + a) = I c (—H a ,7r/2 — a). For this reason, we only present the I c (H a ,a) plots for a in the [7r/2,7r] 
range (with Aa = 15°). The two samples show a quite different a dependence of the normalized pattern width 
A^/A^, as shown in Figs, [sj While the former one is characterized by a weak dependence of Ar(o) with a minimum 
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when the applied field is transverse, the latter one shows a substantial variation with a maximum when the applied 
field is close to be transverse, more specifically, when a ~ 105°. 



C. Annular junctions 

In two recent papers^^, among other things, we have reported on the transverse magnetic diffraction patterns of 
ring-shaped 7V6-based annular JTJs with radii ten times (or more) larger than the Josephson penetration depth. In 
this section we present the results relative to a sample having the mean radius f ~ Aj, namely JJ^G in Table I. Fig(9] 
compares in a combined plot the transverse and the in-plane recorded magnetic diffraction patterns: respectively, the 
open squares referred to top horizontal scale (H z ) and the closed circles referred to the bottom horizontal scale (H x ). 
The vertical logarithmic scale was needed to enhance the plot differences that would be otherwise barely observable 
using a vertical linear scale (a part of the quite different horizontal scales). On a first order of approximation both 
patterns closely follow the zero-order Bessel function behavior; the solid line in Figj9]is the best data fit using Eq.(23) 



with the first critical field as a unique fitting parameter. Fig(9] indicates that for this particular sample a transverse 
field modulates the junction critical current about 1.5 times faster than an in-plane field. In Ref. 17 we have shown 
that this gain increases with the ring diameter and can be even larger than 100. 



As shown in Fig 10 , the annular range width A^ drastically depends on the external field orientation a. In fact, 
although its values for a = and 90° belong to the same order of magnitude, A^ = O.GSA^, we see that A^(a) is 

peaked at olm — 145°, with Aa(o^m) — 100 A^ ~ 160 A^; in other words, when a ~ ckm? the sample is practically 
insensitive to the external magnetic field, and an external field amplitude as large as the irreversible field is required 
to reduce the critical current I c to 67% of its zero field value Iq. 



IV. THEORY 



In order to provide a theoretical interpretation of the experimental data presented in the previous Section, let us 
introduce the spatial normalized units x = x/L and y = y/W with the junction center coinciding with the axis origin. 
Our task is to find out the Josephson phase distribution <p(x : y) over the barrier area (—1 < x < 1 and — 1 < y < 1) 
of a small planar JTJ in a weak oblique applied magnetic field H a . As a preliminary step, before resorting to Eq.Q, 
we need to determine the magnetic field distribution over the barrier area H(x,y). As the Maxwell equations are 
linear in the magnetic field, one can resort to the principle of superposition to calculate the field H. Thus, the effect 
of the oblique field H a can be conveniently split into the sum of the effects of two orthogonal components, that is, 

the in-plane component \/ H% + H% and the transverse one H z . In other words: 



H(x,y)=HHx,y)+H ± (x,y), (7) 

in which H" and are the barrier field distributions induced by in-plane and transverse external fields, respectively. 
As mentioned in the Introduction, it has been traditionally assumed that when the external field lays in the barrier 




FIG. 6: (Color online) Log-log graph of the transverse magnetic pattern I C (H Z ) of the cross square junction. The closed squares 
are the experimental data, while the open circles are the result of computations based on (15). The dashed and solid lines 
are the large-field best fit of the experimental and computed data, respectively, oc H~ 2 and oc H~ x . The inset shows the 
experimental data on linear scales. 



(a) 



(b) 



FIG. 7: (Color online) Tridimensional plots for the recorded oblique magnetic diffraction patterns I c (H a ,a) of two cross- type 
junctions: a) square junction with (3=1 and b) asymmetric junction with (3 = 0.25. The angular separation is Aa = 15°. (For 
these samples I c (H a ,7r/2 + a) = I c (—H ai 7r/2 — a).) 



plane (H z = and H x = 0), it uniformly threads the oxide layer, so that H" = H = H a = {H x , H y ). Today we know 
that this is only true to the first approximation for naked JTJs, since field focusing effects should to be considered 
in planar JTJ structures especially in the case of window junction d 20 ! 21 ! that are surrounded by a passive thick oxide 
layer, the so-called idle regioiW*. Our samples were designed to have the smallest possible idle region, so that field 
focusing effects could be neglected; consequently our theory has been developed under the simplifying assumption 
that the samples are naked. Further, as shown in Figjl] we will only consider magnetic field directions confined to the 
plane specified by the angle a between the applied field and the junction plane [H y = 0). According to Eq. (m), due 
to the linearity of Eq.(fiT), also the phase distribution can be written as the sum of two terms c/)\\(x,y) and (j)±J%,y): 



<l>(x,y) = 4^(x,y) + ^(x,y), 
provided that: V(/>" = ^H" x z and V(j)^ = kH 1 - x z. 

Once y) is known, it will be possible to calculate the junction critical current I c 



(8) 



,14. 



Ioy/ (sin(/>) 2 + (cos< 



(9) 
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FIG. 8: (Color online) Magnetic field range vs. a (in degrees) for the two cross-type junctions of Figs|7| a) (3 = 1 and 
b) (3 = 0.25. The experimental data are presented by closed circles, while the open stars are the result of the calculations 
described in Sec |IV| The insets sketch for each sample its electrode configuration and its orientation with respect to the chosen 
Cartesian coordinates. 
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FIG. 9: (Color online) Comparison of the magnetic patterns recorded for the annular junction quoted in Table I (f « Aj) in a 
transverse (open squares referred to top horizontal scale) and in-plane (closed circles referred to the bottom horizontal scale) 
applied magnetic field. The logarithmic vertical scale helps the data comparison in the lower current range. The solid line 



corresponds to a Bessel-like fit according to Eq. ( 23 ) 



in which the brackets () denote spatial averages over the junction area - 4(/(x, y)) = dx dijf(x, y). 
Being HH = H x St (for naked small JTJs), <t>W is given by Eq.Q that, with our normalization, becomes: 



0" (x, y) = h^ycos a, (10) 

in which ft" = 27r/j J oH a Wd e /<&o, since H a cosa = H x . (If one wants to consider the effect of the idle region, still 
ft" oc H a , but the proportionality constant needs proper correction.) and their corresponding </> x have been found 
in Ref. 9 for naked small JTJs having the most common electrode configurations. Generally speaking, it was found 
(j) 1 - oc H z = H a sin a, so that we are allowed to write: 



<t>^{x,y) = h ± (j) ± (x,y)sina, (11) 

with h 1 - oc H a and (j)±(x,y) containing the spatial part of <j)^{x,y). It is important to stress that the proportionality 
constant between h 1 - and H a is not known a priori, being a non-trivial and still unknown function of several geometrical 




a(degrees) 



FIG. 10: (Color online) Magnetic field range A a vs. a (in degrees) for the annular junction quoted in Table I (f w Aj). The 
experimental data are presented by closed circles, while the op en s tars are the result of calculations described in Sec |IV| The 
solid line arises from the approximate analytical expression Eq.(27). The inset sketches the Lyngby-type annular junction and 
its orientation with respect to the chosen Cartesian coordinates. 
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junction features such as the widths and the thicknesses of the two electrodes, their separation and their configuration^. 
Therefore h 1 - remains the only free parameter when comparing the experimental data to their theoretical counterparts. 
Eq.([8| can be now rewritten in terms of Eqs.(lO) and (11) as: 



</>(x, y) = ft "2/ cos a + h- L (j)j_(x, y) sin a. 

Being both ft" and h ± proportional to the intensity of the applied field H a , their ratio r\ = ft- 1 /ft" depends uniquely 
on geometrical parameters (including the junction magnetic thickness d e ). Therefore, the last equation can be conve- 
niently cast in its final form: 



<t>h,ot(x,y) = h[yc 



- rj(j)±(x, y) sin a] 



(12) 



in which, ft" has simply been renamed ft and the indices ft and a have been added to explicitly indicate that the 
Josephson phase distribution depends on the field strength and orientation. In the remaining part of this Section, we 
will extensively make use of Eq.(12), inserting, for each junction geometrical configuration, the proper expression 
from RefP. 



A. Overlap-type junctions 

From the numerical analysis of the magnetic scalar potential induced in the barrier plane of a planar JTJ by a 
transverse field we were able to derive approximate and simple expressions for the Josephson phase distribution in the 
barrier area that satisfy the Laplace equation (d 2 (j)/dx 2 + d 2 (j)/dy 2 = 0). These heuristic expressions were found to 
be markedly dependent on the junction aspect ratio (3 = L/W. For an overlap-type junction in a unitary transverse 
magnetic field we found: 



<f>±(x,y) « sin y 



cosh (3x 
sinh/3 ' 



(13) 



Inserting Eq.(13) into (12), we observe that 



considering that 



bh, a is an odd function of y, henceforth (sin fa, a ) = 0; furthermore, 
•)h a is an even function of x, the calculation Eq.(|9|) of the magnetic diffraction pattern reduces to: 



I c (h,a)=I dx dy cos fa iOL . (14) 
Jo Jo 

The above integral has been numerically evaluated as a function of the reduced field ft and for several values of a, 
setting Iq = 3.9 mA, (3 = 1 and r] = 0.85 in order to make a comparison with the experimental data for the square 

overlap junction of FigsHl- the calculated I c vs. ft are reported as dotted lines and a proper horizontal scaling was 

II 

chosen to match the in-plane field range = 1.29raT. 

By construction, for a = 0, Eq.(fl4]) returns the Fraunhofer pattern. For the calculated I c (ft), at a quantitative 

level, is only consistent with the experimental data, the discrepancies being more evident for large field values. 
However our calculations can grasp most of the pattern small field features: in particular, they can nicely reproduce 
the dependence of the pattern width A# on a, as far as (3 < 1 (see the open stars in Figs. [5^L-d) . In fact, for (3 = 4, 



to use expression (13) only allows us to reproduce the correct value of olm- 



B. Cross-type junctions 

For a cross- type small naked JTJ with aspect ratio (3, it was heuristically found 9 : 



. _s'mhpx . _sinhy/(3 
^ ^ V)= SmV ^hp + Sm X ^ * (15) 

This approximate expression, beside satisfying the Laplace equation, has the proper symmetry properties required 
by the problem, <p^_{x,y) = cj) l ^{y,x). We will consider first the relevant case of a square cross junction, i.e., (3 = 1. 
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When this is the case, by retaining the first terms in the Taylor expansion of the trigonometric and hyperbolic 
functions, the above expression (15) reduces to: 



</>±(x,y) = xy, 



(16) 



proposed by Miller et a/P^I in 1985. The transverse magnetic pattern of a square cross junction can be computed 
through Eq.([9), by inserting Eq.(15) into Eq.([8| and setting a = n/2. The resulting I C (H Z ) is shown by the open 
squares in Fig 6] Very similar results are obtained, if Eq.(15) is replaced by Eq.(16). For large fields, the calculations 
are well fitted by an inverse proportionality law, I c oc 1/H Z , while the best fit to the experimental data results in 
a quadratically decreasing dependence, I c ex H~ 2 . Again we come to the conclusion that the empirical expression 
taken from RefPis only valid in the small field range, in this case, as far as I C (H Z ) > 0.2Iq. Also for cross junctions 
the I c (H a ,a) were calculated and the a dependencies of A# were extracted and shown by the open stars in Figs{8] 



In evaluating the integral ([9|, the parameter rj in Eq.(12) was chosen to be equal to its experimental counterpart 
A^/A^. As expected, the calculated A#(a) are symmetric with respect to a = 90°. 



C. Annular junctions 

In this section we will examine the static behavior of a small annular JTJ in the presence of an oblique magnetic 
field. Denoting the inner and outer ring radii, respectively, as and r , we assume that the annular junction is 
unidimensional, i.e., the ring mean radius f = (r^ + r Q )/2 is much larger than the ring width Ar = r Q — T{. 
Using polar coordinates r and such that x — rcos0 and y = rsin0, the Josephson magnetic equation ([!]) can be 
split into: 



— = kHq and — = kH v 
or rati 



(17) 



where H r and Hq are the radial and tangential components of the magnetic field in the ring plane, respectively. With 
the annulus unidimensional, we can neglect the radial dependence of the Josephson phase 23 , i.e., <fr(r, 0) = <p(f,0); 
henceforth, 



-KT 



J d6H r (r,( 



(18) 



in which <j>o is an integration constant. By resorting again to the superposition principle, we can readily write the 
analogous of Eqs.([7]) and ((sj) for annular junctions as: 



H r (6) = Hl{6) + H]r{6) 



(19) 



and 



(20) 



provided d^/dO = —nrlll and d^/dO — —nrH^r. 

It is well knowrP^that, when an external field is applied in the plane of an electrically short (r < Aj) annular junction, 
it fully penetrates the barrier, H" = i^x, whose radial component Hr = H x cos6, through (18), leads to: 



(/>H(<9) = /V 1 sin (9, 



(21) 



with = nH x r. As far as 0(0) is an odd (periodic) function the calculation of the maximum critical current reduces 
to the following integration: 



h r 

n Jo 



d6 cos cj>(6). 



(22) 
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Inserting <j> as in (21 ), we obtain the in-plane magnetic modulation pattern: 



(23) 



in which Jo is the zero order Bessel function (of first kind). In deriving (23), it was assumed that the Josephson 
current density is uniform over the barrier area and that no magnetic flux is trapped in between the junction electrodes 
[(/>(# + 27r) = (j)(0)]. The static properties of an annular junction in a transverse field have been numerically investigated 
in Ref. 9 for a Lyngby type annular JTJ obtained by two films having the same widths 24 ; it was found that, in a first 
approximation, H^r sinusoidally depends on resulting in a Bessel-like transverse magnetic pattern. In other words, 
small differences are expected in comparing the shapes of the in-plane and transverse magnetic diffraction patterns of 
an annular junction, as shown by the logarithmic graph in Figj9] However, a small amplitude third ^-harmonic has 
to be added in order to correctly reproduce H^r{6): 



H^r (<9) oc H z (cos (9 + 3(5 cos 3(9) , (24) 
with the coefficient 5 much smaller than unity. The last expression readily provides a more realistic (j)^(0) dependence: 



6 X (<9) = h- L (sin0 + (Ssin30), 



(25) 



with h ± ex H z . Being H x = H a cosa and H z = H a sma, by inserting (21) and (25) into (20), the Josephson phase 



in the presence of an arbitrary oblique magnetic field applied with amplitude H a = H% + if 2 and orientation 
a = aictaii(H x /H z ) becomes: 



i(0) = fe [sin 6 cos a + 7^(sin 6 + S sin 30) sin a] 



(26) 



in which, again , fell h as simply been renamed fe and rj = fe^/fe". In order to reproduce the experimental findings 
reported in Sec III C, we have computed the I c (fe, a), inserting Eq.([26|) in Eq.(22), being still 4>h,a(—6) = — 0h,a(0)- 



The value of r] has been taken from the experimental A^/A\ ratio: tja = 0-63 from Table I, while the value of S 
was determined from the best fit of the experimental transverse magnetic pattern, 5 = 0.02. The results of such 
calculations for several a values are displayed in Fig{lO]with open stars; the agreement with the experimental data 
is excellent. Indeed, the third harmonic correction is mainly needed to reproduce the peak in A^(a); in fact, with 5 
set to zero, the field amplitude dependence in Eq.(26) would be undetermined for a — — arctanl/77, resulting in an 
unphysical independence of I c on H a . However, far enough from this critical angle, we can set 5 = in Eq.(26) and 



the integral (22) trivially results in a Bessel-like behavior with critical field H c or pattern width A ^4 given by: 



A A (a) = 



A 



|cosa + ?7sma| 



(27) 



The last approximate expression, plotted as a the solid line in FigflO} exactly matches the experimental points, 
everywhere except near olm = — arctan( 1/77,4) c± 145°, where it goes to infinity. From measurements not reported 
in this paper, we found that Eq.([27|) can be usefully applied to reproduce also the behavior of long (r > 10Aj) 
unidimensional Lyngby- type annular junctions in the presence of an arbitrary oblique field. 



V. DISCUSSION 



A similar approach can be adopted to describe the behavior of overlap- type junctions having the aspect ratio smaller 
than unity. In fact, being — 1 < x < 1 and — 1 < y < 1, then sin?/ w and, under the assumption /3 < 1, cosh/3x ~ 1. 
Therefore, Eq.(12) simplifies to a linear y dependence: 



(j)(x, y) = fe(cos a + 77 sin a)y, 
resulting in a Fraunhofer-like magnetic diffraction pattern with H c or, with our notation, A#, given by: 
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|cosa + 77 sin a | 



cos a sm a 



x r R 

It's seen that the critical width Ar diverges at a critical angle a c given by: 



(28) 



OLr 



- arctan - 



1 



arctan — + . 



(29) 



In other words, when a = a c , the effect of the transverse component of the applied magnetic field exactly cancels 
that of the in-plane component, so that the junction is virtually insensitive to the magnetic field. In practice, a full 
compensation is never achieved, although A^cxm) can be as large as 100 A^. Eq.(28) is plotted as a solid line in 
Figs{5^i-c and it is seen that it nicely reproduces the experimental data, not only when /3 < 1 (see Figj^), but also 
for (3 = 1 (see Figj5^i). It also closely fits th e results of a long (2W ~ 12 Xj) unidimensional JTJ (see FigJEJ^) for 
which the calculations developed in Sec IV A for electrically small junctions do not apply. Eqs.(27) and ( [28] ) indicate 
that, for small unidimensional annular junctions and small overlap type junctions with small aspect ratio, only two 
measurements are needed to forecast their static behavior in an arbitrary oblique field, that is, the in-plane and the 
transverse magnetic diffraction patterns. 

As far as cross-type samples J Jj^E and JJ#F are concerned, we observe that the experimental data shown, 
respectively, in Fig{5|a) and (b) are affected by a slight asymmetry with respect to a = 90°, in contrast with the 
system symmetry properties (even in the case of electrically long junctions). We explain this symmetry break in 
terms of tiny fabrication misalignments. In our fabrication line, once the base electrode has been etched away, each 
next layer is positioned with an accuracy better than 1 \im (in each direction). This accuracy, although smaller, is 
comparable with the idle-region dimension, resulting in an unavoidable, albeit small, imperfection in positioning the 
junction area exactly in the centers of the films. Since the screening currents mainly flow along the electrode borders, 
the misalignment effect increases with the junction side dimension. This explains why the observed asymmetry is 
larger for junction #F. 

For the Lyngby-type annular junction our theory reproduces the experimental data in a more than satisfactory 
fashion. However, the situation might be not so good for annular junctions made by electrodes of unequal widths as, 
for example, those used in RefP^, which require also the introduction of the second ^-harmonic in Eq.(24). However, 
when (j)^{9) is as in (25), exploiting the trigonometric equivalence sin3x = sinx(l + 2 cos 2x), the transverse magnetic 
pattern can be analytically shown to be given by the following even expression: 



I c (h)\) = / |J (/i ll ) + 2(5J 2 (/ill)|, 



in which J2 is the second-order Bessel function and 5 was assumed to be much smaller than unity. Similarly, I c (h) 
can be analytically worked out in the more general case of Eq.(26), as far as a is far from %, more precisely, when 
S « I + I/77A tan a. 

The weak point of our theoretical approach based on the superposition principle is the indetermination of the 
parameter h 1 - introduced in Eq.(ll). The knowledge of the magnetic field actually introduced into the barrier for a 



given transverse field H z requires a careful experimental investigation with samples having a given junction geometry 
and different geometrical parameter of the connecting electrodes. However, while the superconducting film widths 
can be easily varied, it is difficult to realize samples with much different film thicknesses. Consequentially, when the 
shape of the I C (H Z ) is known a priori, the effective normalized field h 1 - remains to be determined from the direct 
measure of the transverse magnetic diffraction pattern. 



VI. CONCLUSION 



In this paper we examined the static properties of small planar Josephson tunnel junction in presence of a uniform 
external field applied at an arbitrary angle with respect to the barrier plane. This topic has been considered from 
both the experimental and theoretical point of view. We have presented the recorded oblique magnetic diffraction 
patterns of junctions with the most common electrode configuration, namely overlap, cross and annular geometries. 
These data, beside being original by themselves, also served as a test for the theoretical analysis of small JTJs in a 
purely transverse field that we recently proposed 9 . Further, by invoking the superposition principle, the findings of 
RefP in a transverse field were combined with the classical knowledge for a JTJ in a parallel field to provide a general 
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theory for any arbitrary oblique field. We stress that the theory has been developed assuming that junction was 
electrically small, naked, geometrically perfect and made with purely diamagnetic superconducting films (Al = and 
H c = oo). Although our samples satisfied these conditions only to a rough approximation, the agreement between 
the experimental and theoretical results is more than satisfactory especially for annular junctions. For rectangular 
samples the theoretical approach is again very good only for small values of the applied field and of the junction aspect 
ratio. In the other cases, our theoretical predictions fail to provide a correct description. This was to be expected 
because the approximate analytical expressions heuristically found in RefP were already observed to have the largest 
relative error near the junction corners. As already reported elsewhere^, the importance of the Josephson phase at 
the junction corners grows with the amplitude of the external field. 

The main message of this paper is that even a small transverse field (which has been largely ignored in the past) 
can strongly influence the magnetic interference patterns. We explore the implications of this result in supposing 
systematic errors in previous experiments and in proposing new possible applications. Since in most of the applications 
the external field needed to modulate the critical current of a planar JTJ is applied in the barrier plane by means of a 
long solenoid or Helmholtz coil pairs, we want to stress the importance of the alignment of the solenoid (or coil) axis 
with the junction plane. It was believed that any possible tiny angular misalignment 5 a between the coil axis and 
the junction plane would result in only second-order errors, being H x = (1 — 0.b5a 2 )H a . However, as we have shown, 
this is not necessarily true, since the unwanted (and often unknown) small transverse field component H z = 5aH a 
might have an effect comparable or even larger than that of the in-plane component, if the junction critical angle a c 
is close to 180°. For example, in the case of sample JJ#C of Table I, a misalignment 5a c± —4° would result in a 
systematic error of more than one order of magnitude for the measurement of the first critical field. The consequence 
of the coil misalignment might have been underestimated, if not ignored, in many previous experiments dealing with 
JTJs in an external magnetic field, including those in which the exact knowledge of the field in the barrier plane is 
of capital importance. As a corollary, it also follows that in shielding a cryoprobe the same care has to be taken to 
minimize both the in-plane and the transverse stray fields. 

Furthermore, in planar SQUID applications the magnetic field to be measured is applied perpendicular to the SQUID 
loop and its effect on the junction critical currents has never been considered. However, it might not be negligible, 
especially when the field is large and the junction(s) is placed close to the borders of the superconducting electrodes 
where the induced screening currents are larger - this is the case of step-edge or ramp- type junction J 25 | 26 l 

In this paper we have shown that for a given junction geometry the response to an externally applied magnetic field 
drastically depends on the field orientation, a property that might be exploited to design angle resolving instruments. 
Furthermore, considering that this response is different for JTJs having different geometries it makes possible to 
design multijunction chips in which, for a given applied magnetic field, the critical current of some junctions is almost 
completely suppressed while that of other junctions remains unaffected. Our findings also suggest that two (or more) 
independent magnetic fields with different amplitudes and orientations can be applied to multijunction chips in order 
to obtained the proper critical current suppr essio n required for samples having different geometrical configurations. 
By using Eq.Q, the theory developed in Sec IV for a uniform oblique field applied in the y = plane, can be easily 
extended to the most general case in which all three field components are non-zero. In other words, the oblique 
magnetic diffraction pattern can be theoretically predicted for the most common junction configurations, as far as the 
junction dimensions are smaller than Josephson penetration length. The case of long JTJ in a transverse field still 
remains an open question since it requires the solution of a tridimensional magnetostatic problem in the presence of 
external (non-linear) currents 2 ^. 
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Note 



Soon after this paper was accepted for publication in Phys. Rev. B we become aware of a paper by Heinsohn et 
al. entitled "Effect of the magnetic-field orientation on the modulation period of the critical current of ramp- type 
Josephson junctions " (see Ref.27). The authors report several experimental magnetic diffraction patterns of both low 
and high temperature ramp- type Josephson tunnel junctions. As first found by Rosestein and Chen, they observed 
that the pattern modulation period drastically changes with the external field orientation. 
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